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Consider  a  two  dimensional,  isotropic,  canonical  minimum  scattering* 
antenna  in  a  homogeneous,  lossless  medium  set  in  the  plane  y  ■  0.  This 
antenna  produces  a  two  dimensional  E-type  field  i.e. ,  E  parallel  to  the 
z-axis.  (see  figure  1.)  Taking  E(x,y)  ■  E(r,+)  where  x  “  rcos  + 
and  y  “  r  sin  $  ,  the  spatial  Fourier  representation  is 


E(x,y)  -  E(x,y)z  -  f  |(u)e“lll*"lvy< 

~  J  * 

/'  A  -ik*r 

E(u)e  '  ”du  i 

* 


where  k»ux  +  v^,t«k/|k|,  v  » 


u2  |u |  «  6 

Li/2-V  W  >  e 


with 


B  ■  2n/X  and  implicit  time  dependence#  e+iwt ,  assumed  and  suppressed. 
Note  the  choice  of  cuts  for  v  is  consistent  with  the  restriction 
y  >  0  for  outwards  traveling  waves.  On  setting  u  ■  8  cos  (, 
v  ■  B  sin  E  transforms  as 


E(x,y)  -  E(r,4)  •  B  I  EW"*6*  •  sin  i  d( 

•'s 

where  S  is  a  Soamerfeld  contour  in  the  (-plane  -  see  figure  2.  Deforming 
S  to  a  steepest  descent  path  in  the  region  of  convergence  for  E  gives 
the  asymptotic 


Ve  use  freely  the  previous  study  report  [1}  by  Professor  Prancea- 
chsttl.  Constraints  ae  "invisibility",  i.e.,  antennae  that  do  not  scat¬ 
ter  with  particular  loadings  of  ports  -  "Canonical  Minimum  Scattering" 
being  as  such  with  ports  open  circuited  -  are  discussed  at  greet  length 
therein. 


8.D.  Path 

Figure  2  Sommerfeld  Path 
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E(r  »  1,4)  -  -0  sin  4  1(6  cos  4)  a"4®' 

It  is  possible  to  define  a  noraslixed  antenna  curr«t(  I f  such  as  to 
require  (1)  reciprocity,  f(4)  ■  *r#c(4),  where  i»  the  trans¬ 

mitting  (receiving)  effective  height  per  unit  length.  This  inposes  the 
following  constraint  on  spectrum  and  pattern  functions  for  real  valuea 
of  the  observation  angle  4  €  10, r] 

E(6  cos  4)  - 

4w 

The  Identity  theorem  provides  for  the  coincidence  of  the  above  quantities 
throughout  the  coonon  domains  of  analyticity.  Thus, as  f(u)  -f ( 0cor  O, 
we  know  the  spectrum  |(u)  if  we  can  find  the  analytic  continuation  of 
the  pattern  function,  f(u),  from  the  visible  region  u  4  [-6,8]  through 
all  regions  of  coamon  analyticity.  As  f(u/  6  8(u)  are  often  entire 
functions,  this  region  le  generally  all  of  C.  Lastly, we  note  that  the 
Fourier  transform  of  the  magnetic  field  asaociatad  with  this  E-polarlxed 
field,  fi(u)  has  the  following  properties} 

1)  H(x,y)  -  H(x,y)s  •  I  B(u)s  e"ik*rdu 
*  ■'* 

2)  c  g(«)  •  *  *  |(«) 


With  these  resulte  it  is  eaey  to  find  the  total  power,  Pfcot 
into  the  region  y  >,  0  by  the  antenna 


I(*,0)  » 

CB'x,0)  - 


(t  x  s>E(u)e*’1*,*du 


injected 
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ptot-i  f  *  * »’ *  *■  •  $  w  /  *<•>»*<«>.*. 

"'ll  ■'ll 

.» [  i!si£4ifc 

5"ii?u  4  * 

gm  9 

Frequently  ?tot  »  ^  |l|  where  Z  -  R  +  iX  ,  the  entenne  system  input 
lmpedance(e,g.  "one  mode  slots"  In  ground  plane  or  infinite  wires)* 

Here  one  has  a  relation  between  a  measurable  quantity  2(u)  and  f(u) 
which  we  shall  exploit  presently.  Rote  this  constraint  applies  to 
both  "visible"  and  "invisible"  values  of  f(u). 

Finally  we  give  a  formal  expression  for  the  mutual  coupling  between 
antennae.  (See  Figure  3.)  Reference  lines  of  Antennae:  1  -  (0.0). 

II  “  (xo,yo)  -  (r0,dQ).  Take  the  plane-wave  spectrum  of  antenna  I  with 
antenna  II  open  circuited  faote  invisibility  assumption  for  II) »  then 
the  spectrum  of  the  radiated  field  is: 

*!<«>  *  ^  * j(C)di . 


The  open  circuit  voltage  at  antenna  XI  due  to  the  spectral  component 
Ej.(u)du  is 

-lux  -ivy 

a 

dV21  * - - -  du  fI(*)fII(T) 

A*v 

A  ^  a 

Now  for  given  C,  k  •  cosfcx  +el*  \  y  whereby  the  radiation  Impinges  on 
II  at  angle  {  ♦  i  (see  Figure  3.)  thus  by  superposition: 


BC  1,  /*  fx(«fIX(C  ♦  w)  - 

hi *  —  J - ; - « 


•i™.  -  lvv 


Figure  3  Mutual  Coupling  Geo 
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V_.  mutual  Impedance  . 
Z_.  »  ~ —  -  between  2  C.M.S.  2a 
X1  antennae. 


fI(0fIIU  *  "ie'^V^o 
v 


du 


-/ 


txtt)f2u 


-18r  coa(d  -O 
+  ir)a  0  °  dC 


by  the  pxevioua  change  of  variable. 

Now  we  assume  the  antennae  to  be  identical  (e.g.  In  an  array  of 
slot  radiators  etc.)  whence  these  expressions  simplify  for  real  excita¬ 
tions: 


*<“>  •  £  /  note  -u  *  B  cos(C  +  w) 


If  u  ■  6  cos  €  then  E(-u)  ■  !(u)*  where  E(x(0)  is  the  "antenna 
current"  (equivalent)  excitation  at  y  »  0.  Therefore 


E2X -if  /  - 

y,  v 


d«-4l 


-lux 

P  P(u)e  0 
&  A2  -u2 


du 


where  P(u)  Is  the  conventional  power  density  with  u  ■  I  cos  C  and 
the  antenna  system,  again,  constrained  to  y  ■  0.  If  we  force  the 
sources  to  be  of  finite  extent  the  pattern  functions  become  band-limited 
l.e. 

f(u)  -  f (0  eos  O  -  /  ei(a  °.xJ(x)dx 

*'-e 

where  the  excitation,  j(x).  Is  nonsero  x  <  [-a, a]  (J  -  J(x)6(y)t  la 
an  equivalent  current.)  The  sequel  considers  the  utility  of 

this  constraint  of  band -limitation  and  the  possibility  of  relaxing  it 
as  well  as  the  requirement  of  identicellity  of  radiators  previously  Im¬ 
posed.  We  should  state  here  that  we  are  attempting  to  calculate  the 
mutual  couplings  between  aerials  by  means  of  the  previously  deduced 
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Integrals  Which  require  knowledge  of  the  pattern  or  power  pattern 
function  over  Loth  visible  and  invisible  values  of  argument.  At  hand 
are  the  values  of  the  power  and  pattern  function  in  the  visible 
[-0  £  u  £  0]  region  and  possibly  the  input  impedance  over  soeie  range 
of  frequencies*  The  required  nuaerical  analytic  continuation  to  invis¬ 
ible  values  is  a  very  delicate  process  (a  point  to  be  amplified  below) 
and  requires  much  constraint  to  provide  practicability.  In  that  one 
nay  approximate  any  entire  function  by  a  sequence  of  band -United 
functions*  the  constraint  of  band -Halted  proxlaates  to  a  given  pattern 
function  is  not  too  onerous.  While  the  general  problems  of  approxima¬ 
tion  in  this  context  will  be  considered  extensively  below,  the  virtue 
of  this  scheme  will  be  displayed  now. 

Consider  the  power  like  quantity  fg(u)f 22<-«)*  Mote  that  given 
*1**2  *  *  Ti>*  ProducC  quantity,  f  j(u)f  jj(-u),  la  also  band-limited  as 
follows  from  the  convolution  theorem*.  In  fact,  f  j(u)f IX(-u)  C  B^  » 

In  the  case  of  identical  aerials  this  f(«)fM*  P(u)  -  the  convention¬ 
al  power  pattern  function, which  may  be  measured  with  eons  accuracy  there¬ 
by  constituting  a  reasonable  candidate  for  numerical  analytic  continua¬ 
tion.  (The  inaccuracy  of  phase  naasuriant  would  seem  to  forbid  an 
attempt  at  continuation  of  the  pattern  function  itself.)  Let  D  be  the 
distance  between  aerials*  reference  lines  than 


Space  of  band-limited  functions  l.e.  f(t)  •  J  e*"  g(w)du  where 
g(»)  la  an  t)  function  i.e.  "square  lategrable." 

*F1(»)P2(»)  *■*’  fl(t)*f2(t)  where  f(t)  ■ /Re4“tF(w)|y  etc. 
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Then  P(u)  *  £(u)»f  ( -c)  •  jf(u)J 


2a 


2  *  (g(x)*g(-x)](x)  -  G(x) 
x  £  [-2a,  2a] . 

■  0  otherwise 


Hence  P(u)  ■  J  G(x)  e“luxdx.  Expending  G(x)  in  an  exponential  series 

2a 


-2a 

in  I -2a,  2a];  G(x)  •  IC  eln1^c/2•  with  C 

v  n  n 


-  £  f  G(x)  e^^dx 


-2a 


P  (^~)  ■  So  one  gets  the  "saapllng”  theorea  representation: 
2a  2a 

«»)  -  f i  .  i  c  f  -  »> 

H  J  ®  J  “ 

-2a  .2* 


-2a 

.  t  Sin(2au-ns) 

J,  *'2a'  2au-n* 


How  insert  this  result  into  the  coupling  integral  and  invoke  the  "con¬ 
volution  theorea" 


Lt=t  du  ■  v®  I 

”"/g2-u2  *  J  2 


4w" 


du 


£  P(SJ)l(u,B,D) 

at  •• 


Recall  the  following  fourier  transfora  pairs 


1)  Sln(2au-nw) 
2au-nw 


.ixn*/2a 


xc  [-2a,  2a] 


■  0  otherwise 


2)  and  for  our  choice  of  cuts  (ef.  p.  1). 


hp^ 


j 

! 

c 

» ! 

f  *  ■ 

) 


hJ2)  (e|x|) 


/(*2-u2 


hence 


2a 

(n.B.D)  -  gj  J[lnxx/2e  R(2)  tgjD  +  x|j  dx>  Now  recall  further 


-2e 


(2) 


the  Integral  representation  for  H'  over  finite  limits 


.2  t  j 


2r  Cot  6 


d0  (Re  a  >  0} 


Sin0  /cos 6 


This  may  be  used  In  the  expression  for  I(n,0,D)  and  limits  of  inte- 


Lti 


gration  exchanged  to  reduce  161,8,0)  to  a  concatenation  of  elementary 
functions  readily  amenable  to  evaluation  with  1)  fortran  complex 
algebra  2)  any  simple  numerical  integration  scheme,  e.g. 


for  D  >  2a  ve  have 
*/2 


I(n,B,D)-  2ie 


•180  f  s-i9/2  e 

J  ■  -f  n  A 


-2BD  Cot  8  sin(ns-2aB44gla»0ote) 


de 


sin  6  /cos  0 


(nir-2afi+41aBCot  0) 


with  a  somewhat  more  tedious  form  for  0  <  2a. 

So  for  reciprocal  antennae  with  real  excitation  the  first  order  (exact 
C.M.S.)  mutual  coupling  may  be  expressed 


J  /w 


e'luI>du  -  2>(g)  !<«.•.»>  »ote  that 


-  B 


In  the  integral  formulation  the  reactive  part  of  the  coupling  is  ex¬ 
tracted  from  the  Integration  over  ehe  invisible  (8  <  |uf)  range  of  the 
pattern.  Likewise  the  radiative  portion  of  the  coupling  from  the 
visible  portion,  to  see  this  note 
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8 


‘  *  /. 


P(u)« 


-iuD 


-iuD 


/  a«ul_d» 

l" I*8  JF77~ 


where  P(u)  is  In  general  an  even  function  hence  the  firat  integral  ia 

obviously  pure  real  aa  Is  the  aecond.  However,  we  have  reduced  this 

quadrature  to  a  series  whose  nth  term  has  both  real  and  imaginary  parts. 

Thus  the  "sampling  transformation"  displays  the  contribution  of  both 

visible  and  invisible  parts  of  the  pattern  to  the  total  mutual  impedance. 

It  should  be  pointed  out  that  P(^)  end  I(n,B,D)  are  both  "Fourier- 

serles"  coefficients  of  G(x)  and  +  x  |),  hence  decay  as 

O(^)  and  give  a  somewhat  tardy  convergence  of  0(— — )  to  the  mutual 

n 

impedance  series.  It  is  possible  to  Improve  this  result  somewhat  as  we 

shall  show  later.  Now  however,  the  problem  at  hand! 

How  can  one  extract  the  first  few  P(^i)  from  the  know  P(0) 

and  P(u)  for  u  C  [-8,8]  and  thereby 'estimate  the  coupling  for  two 

identical  aerials?  How  may  one  analytically  continue  from  the  segment 

the  band  United  function  P(u)?  He  should  point  out  that,  depending  on  the 

—2a  2a 

antenna  structure(vislble  range  -  dimensions  of  antenna  *  (— y  *  -y)) 

P(0)  P(-|“)  may  be  know.  In  fact  for  highly  directive  antennae,  in 

which  case  the  coupling  question  is  somewhat  moot,  the  extensive  nature  of 

+  nw 

the  visible  range  contains  sufficient  sample  points  (-8  <  -jj-  <  B)  to 
yield  a  vary  accurate  picture  of  the  coupling.  Exclusive  thereof,  we  are 
faced  with  the  solution  of  the  so-called  "ill-posed  problem"  of  "numerical 
analytic  continuation"  of  band-limited  functions.  The  following  discussion 
considers  numerous  naive  but  tractable  solutions  to  this  problem  and  at 
length  some  sophisticated  yet  awkward  methods  of  solution.  Ve  shall 
require  also  some  information  of  the  sampling  theorem  and  its  cognates 


*.  -•  ii"r‘ii’vii  -I •  rr*-n-  I  i  in"  r  ir  ""f ;'?  *  V  sabWBUfjHMfcafe  A  »i  nr  f  liTi  f 
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to  speed  the  convergences  of  the  various  series  Involved.  We  commence 
with  a  general  discussion  of  numerical  analytic  continuation. 

Analytic  continuation,  cornerstone  of  complex  analysis,  is  not  directly 
numerical  in  character.  The  usual  context  is  that  of  a  "magic  wand"  waved 
over  a  given  functional  representation  known  valid  In  a  restricted  region 
but  convergent  in  a  larger  region.  The  agreement  of  these  representa¬ 
tions  on  any  dense  set  (just  a  limit  point!)  provides  for  total  agree- 
ment  in  any  coanon  domain  of  analyticity.  Hence  one"analytically  contin¬ 
ues",  the  function  from  some  small  set,  with  reckless  abandon,  throughout 
C  or  so  it  seems.  In  the  case  of  an  exact  agreement.  e.g.  two  integral 
representations  of  T(z)  ,  analytic  continuation  is  a  deft  tool.  Numeri¬ 
cal  (experimental)  attempts,  however,  are  characterised  by  the  following 
playful  Theorems  [2],  used  by  Atkinson,  quoted  below. 

Divine  Theorem:  If  one  knows  an  analytic  function  on  a  segment  of 
a  line.  Inside  its  domain  of  analyticity,  one  knows  it  throughout  this 
*  domain  of  analyticity. 


Diabolical  Theorem:  If  the  function  is  not  known  exactly  on  the 
line  segment,  but  only  within  an  error  corridor  of  width  2e  ,  then  the 
uncertainty  in  the  continued  function  is  such  that  its  value  at  any  given 
point,  in  the  domain  of  analyticity,  can  be  any  number  whatsoever,  and 
this  for  any  e  ,  no  matter  how  snail. 

All  seams,  but  is  not,  lost.  It  turns  out  that  "global  bounds"  on 
the  function  to  be  continued  [3]  will  restore  the  "continuous  data  dep¬ 
endence"  lost  in  the  process.  However  the  restrictions  suitable  to  a 
given  problem  are  ad  hoc.  The  sequel  will  examine  two  "almost-bast-possible 
methods  of  analytic  continuation  based  on  the  Miller -Tikhonov  regular is at ion 
schemes.  This  exposition  is  borrowed  with  modification  from  recent  optics 


-  ■  *  -  -  *r  i  -|  *•  ntgttdirtto’lfcr 
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literature  of  Inage  extrapolation  and  object  restoration*  These  method s, 
which  utilize  the  eigenfunctions  of  the  finite  Fourier  transform,  the 
Prolate-Spheroidal  wave  functions*  are  quite  cumbersome  to  employ.  They 
may  be  used,  eventually,  as  benchaarks  to  compare  with  proposed  "sampling" 
proximate*  to  be  introduced  shortly.  It  should  be  remarked  that  we  are 
here  concerned  with  the  analytic  continuation  of  band-limited  functions 
l.e.  given  f(u)  -  f  eiuxg(x)dx  and  f(u),  u  €  [-0,6]  find  f(u) 
for  u  i  [-0,6];  or  (harder)  find  g(x)Vx  €  [-a,a]  in  the  above  Fredholm 
equation  of  the  first  kind,  whence  one  may  compute  f(u)  from  the  trans¬ 
form  relationship  above.  It  is  this  duality  between  the  analytic  continu¬ 
ation  of  band-limited  functions  and  the  Fredholm  equation  that  partially 
motivates  our  requirement  of  band-limitation.  The  general  problem  of 
analytic  continuation,  by  means  of  the  Cauchy  theorem,  can  be  expressed 
as  the  solution  of  a  singular  integral  equation  [3]  which  is  an  insidious 
numeric  task,  whereas,  the  constraint  of  band-limitation  requires  the 
"solution"  tj  a  regular  integral  equation  of  the  first  kind  -  a  simpler 
task  well  considered  in  the  literature.  Let  us  consider  a  band-limited 
function  f(u)  with  f(u )•/  eluxF(x)dx  «  B  ,  .  First  note  [4]  the 


sequence 


eln(2au-uir) 


2au  -  ns 


-2a 

is  a  complete,  orthogonal  sequence  in  X  all 


of  whose  terms  are  band-limited,  C  ^  .  If  one  were  to  consider  the 
distribution  due  to  various  currents  at  y  •  0  in  the  previous  geometry 
with  the  insistence  of  piecewise  constant  currents  (thus  approximating 
a  given  smooth  distribution)  the  pattern  would  be  Of  thla  form.  Bacall 
we  have  previously  the  bend-limited  power  pattern  ae  a  sampling  series 


*»  - 1  »© 


We  might  ask  the  following : 
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1)  Is  It  possible  to  spproxlaate  P(u)  In  the  visible  renge  with 

a  "sampling  series"  of  unknown  coefficients  with  the  constraint  of  antenna 
self-impedance  Batched  as  measured? 

2)  Is  it  possible  to  solve  for  P(^j)  ••  •  truncated)  set  of 

linear  equations  that  we  implicitly  defina  through  the  data  as  follows: 

A.  Record  P(Uj)  for  l-B,Bl  for  all  1. 


B.  Solve  the  truncated  set  of  equations 


N  .  sln(2au .  -  jir) 

*.  Kut)  P(£)  —  i  jV  1  -  1.  « 


for  P(^j-),  j  -  -N,  N.  (recall  P(0)  know) 


C.  Impose  the  impedance  constraint  on  these  near  singular  linear 
equations. 


3)  Is  it  possible  to  hasten  the  convergence  of  the  "sampling  aeries' 


and  provide  for  solution  of  fewer  coefficients,  thus  taking  advantage  of 
our  ability  to  sample  the  visible  interval  at  will?  (e.g.  find  a  deri¬ 
vative  of  P(u)) 

All  questions  are  answered  in  the  affirmative,  and  involve  novel 
techniques  that  circumvent  the  previously  required  usage  of  the  Prolate 
Spheroidal  Harmonics.  We  shall  investigate  each  in  turn. 


1)  First,  one  requires  a  criterion  of  fidelity.  Without  loss  of 

generality,  we  take  a  mean-square  (i2(-2a,2a)>  setting.  Taka 

||  * 

.  sln(2au  -  aw) 

^  mi  (2au  -  nr)  ,  with  ao  ■  P(0)  known,  as  the  ssmpling 

XB  -  2 

(F(u)  -  P(u))  du  ■  minimum,  subjsct  to  tbs 

l  (Sj) 
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constraint  2 ^  is  given.  Recall  2^  stay  be  expresaad  la  terns  of 
the  pattern  function  as  previously  shown.  To  wit 


2la  -  R  +  !X 


■n  r  -a 


B2  -  u 


m  f  Jsal—  du 


The  quadratures  nay  be  found  by  direct  integration  of  the  above 


quantities  or  the  previous  convolution  approach  i.e. 


zla  "  R  +  « 


"  If  X 


* 


du  "  If 


.sln(2u -nirh  du 
'  2au  -  n*  '  - - 


r^u2 


■  #  £„  «.2<«.».o) 

n»-« 


where 


Kn.B.O)  -  fU  £  «“ln1rx/2*  H<2)(B|x|)dx 


(even  4  odd) 


(even) 


which,  again,  nay  be  reduced  to  a  ld  integration  of  elenentary 

functions  by  the  prior  Integral  representation  for  the  Henkel  function. 

N 

Ne  have  two  constraints  on  (a.)  „  : 

-  n  n«-N 

"’If  f  •n*«<l(ii,e.O»  X-$  f  .nto(I(n,B.O» 

H"“n 


It  is  now  a  slaple  natter  to  use  Lagrange  nultlpliers  to  solve  for 

(a.)  such  that  /  (P(u)  -  P(u))2du  ■  nlnlaun  given  X,R.  The 

1  •'-B  (a.) 


functional  of  {a^}  to  be  nlnlaised: 


JU‘’  ‘  £  %  •»  -  '<«»’*■  *  ‘i  §  v; + *2  5  vs 
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•  . 


where  I' 
n 

!fx-x. 

require* 


Re  I(n,0,O)  etc.  and  •»  R  »  -  £  a  1*  , 

N  3J{a.)  C6  °  tS  “  n 

Ta  I"  .  t - —  -  0  for  all  a4,  i  -  +  1,  +  2 . +  M 

-N  n  n  3*i  1  ~  ~ 


N 

2  £ 

n»-N 


ain  2au  -  nn 
2au  -  nir 


dv  -  2  p(0)  du 

(2au  -  uir)  av  J_&  1 W  <2au-iir) 


+  XI*  +  x2i”  -  0 


The  last  2N  +  2  equations  may  now  be  aolvad  for  the  2N  +  2  unknowns 
{aj}^_  n,  X^,  X2»  To  minimize  sampling  of  P(u)  in  evaluating  the 
second  integral  one  could  use,  say,  a  Gauss-Legendre  quadrature.  Note 
that  the  finite  range  and  analytic  nature  of  the  integrands  makes  their 
evaluation  a  simple  and  Inexpensive  task.  Results  numeric  of  this 
procedure  are  contained  in  the  first  appendix. 


2)  The  possibility  of  using  the  continuum  of  pattern  values  in  the 
visible  to  extract  a  continuation  through  the  sampling  theorem  is  not  new 
[5].  One  may  measure  at  will  and  solve  the  set  of  linear  equations  des¬ 
cribed  above.  Goodman's  classic  text  on  optics  [6]  considers  this  proc¬ 
edure  and  provides  it  blessings.  Zn  point  of  fact,  one  is  apt  to  get 
order  of  magnitude  results  or  better  In  its  blind  application.  However, 
the  linear  set  is  not,  to  begin  with,  favorably  conditioned  and  any 
attempt  to  continue  "very  far"  will  lead  to  trouble.  The  problem  of 
constraints  or  regularization  in  this  connection  has  not  been  explored  in 
the  literature,  nor  has  the  very  specialised  nature  of  the  linear  equa¬ 
tions  involved  (to  wit  they  have  the  form  of  the  "double-alternant").  He 
Intend  to  pursue  these  toplca  at  length,  and  while  results  to  date  are 
incomplete,  we  give  cause  to  further  these  activities.  Part  of  the 
problem  stems  from  the  slowly  convergent  nature  of  the  sampling  series. 
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This  we  shall  examine  in  some  detail,  offering  "quicker"  sampling 
proxlaates  as  naans  to  the  required  coefficients.  Lastly,  we  indicate 
how  the  transformed  expressions  for  nutual  coupling  might  be  improved 
through  this  device. 

Recall  that  we  wish  to  solve 


M  ,  sin(2au.  -  jw) 

»<V  •  E  Kg) - 4 —  , 

*  j— H  2au1-jir 


1  -  1,2H 


for  P(^)  J  *  "N,...,||,J  f*  0  (which  is  known,  P(0>)  by  accurate 

samplings  of  P(u.)  for  u,  [-0,6]  subject  to  the  self -impedance 

N  N 

constraints,  truncated  as  X_  -  J2  P(5^)IW  ,  R  ■  12  *(5^)1'  with 

0  -N  2*  “  °  2a  n 

I"  etc.  as  previously  defined.  For  the  moment  we  ignore  these  constr¬ 
aints  and  consider  just  the  set  of  linear  equations  above.  As  previously 
remarked,  direct  solution  thereof  produces  "trash"  for  N  of  any  alxe. 
Out  of  curiosity,  the  least-squares  (normal)  solution  was  also  considered 
for  a  few  test  functions  (Appendix  B,  pt.  1)  l.e. 


If  A  5  •  where  b  «— *■  (PO^)};  5  «— ►  (K^))  £  ►  (- 


slnflauj-j  w) 
2au4- J  * 


) 


T  T 

the  nonsingular  least  squares  problem  A  A  x'  •  A  b  has  x’  as  solution 
to  the  problem  ||  Ax'  “  b||2-  •  minimum,  uniquely.  The  results  (Appendix  B) 
indicate  order  of  magnitude  agreement  "in  trend"  with  the  true  coefficients. 
(This  will  be  seen  consistent  with  Tikhonov  regular isat ion  of  the  lllposed 
linear  set.)  Furthermore  these  truncated  equations  are  of  a  peculiar  form. 
Previous  investigators  (7,8]  have  considered  similar  sets  and  noticed 
closed  form  solutions  may  be  readily  obtained  from  Cramer's  rule  as  such 
equations  have  the  form  of  the  double-alternant  [9]. 

To  see  this  we  renormal lx e*  the  aquations  as  follows: 


assume  F(0)  unknown  to  simplify  indices. 
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sln(2aui  -  Jt)  -  [tin  laUjH-l)**  lat  P(^j)  *  (-l)^P(^y)  -  lj  , 
P(ut) 

8 ij>  2au  *  Bi  »  Yj  •  I(j  -  1)  -  HJir,  =  Pau^.  Thua  we  obtain 


B 


i 


2N+1 

L 

J-i 


i  -  1,2K  +  1 
j  •  1,2N  +  1 


Now  the  determinant  of  C  may  be  found  In  simple  closed  form.  Some 
more  notation  [9]  la  required  however. 

Define  the  "difference-product"  of  the  array 

as  follows 


S1,2(a1#e2,...a  )  •  T^f  TT  (*.  -  s.) 

12  n  i+2  J-l  1  J 


then  (9] 


I  Cl  -  Det | C | 


2N+1 


A.tr 


2N+1 

TT  «K 

K-l  * 


where  is  a  cross  product 
rule  on  (*  C  A  and  taking 


IS(|< 


D_  ■  TT  (*r  *  UslnB  Cramer’s 

a  i-1  R  1 

as  the  deletion  of  the  ith  column. 


S 
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jth  row  of  | C |  In  the  usual  expansion  about  minors,  gives  the  following 
explicit  form  for  the  solution  of  the  truncated  (infinite)  set: 

!-$i-  *t  ’  $  w,3+1  ^  -*J 

where  are  the  (measured)  r.h.s.  of  the  normalised  equation  set  allow¬ 
ing  one  thereby 

1)  to  see  how  an  error  In  (Bj)  might.  pt%g* «#,?*• 

2)  the  explicit  calculation  of  the  pi>so*;'  ^iutlon  and  suggestive 
perhaps  of  a  new  analytic-product  form  £ep$'«fwt»$fitlve  of  the  band-limited 
function. 

3)  to  settle  the  question  of  what,  if  any,  constraints  on  (B^J 
are  necessary  to  obtain  the  solution  of  these  equetlons  [e.g.  5]  in  the 
presence  of  noise  or  otherwise  (Linfoot  and  Shephard  [8]  found  In  a  simi¬ 
lar  but  not  applicable  set  explicit  constraints  on  the  various 

(Bj),  { 6^  -  «tc.  dictating  convergence.)  We  are 

presently  attempting  to  resolve  these  questions  and  can  provide  here  only 
a  few  Initial  observations. 

As  the  row-column  deletions  leave  yet  another  alternant  form 

lit,  I 

- •  A. .  is  mildly  tractable,  just  delete  y.  and  6.  from  the  ex- 

|c|  X  3 

presslon  for  |C|  and  divide  obtaining 


L  m 

ji 


iB+l  ZNtI 

_ ~  v  U  (gi  -  v  * 1 _ 

i-l  2Wfl 

TrZN(i-l)!(2W  +  l-i)!  («.  -y.)  TT  <«4  "  V  TT  <«r  “  «*> 

3  1  K-l  3  *  KM+1  *  3 
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hence 


2N+1 

.  U  (<5K  "  Yi)(“1)  2Si1  I-I  (gj  ~  VB1 _ 

1  "  2N  i-1  44-  2*4 

7rZN<l-l)!(2H  +  l-l)!  31  («.  -  Y.)  TT  («4  “  «J  FT  <«r  -  «*> 

3  1  K-l  3  *  r-i+i  *  3 


2N+ 1 

L 


It  le  now  neceeeery  to  assuae  some  regularity  on  the  sample  points  <  . 

& 


“  si+i  -  5i 


A,  for  all  i,  the  monster  simplifies  thus 


m1  rt«k  -  v  ^  <-)3>1  sPf<«1  -  vk> 

^  _  k-l  k  1  £  _ 3  k»l  3 

1  A2Nir2N(l-l)f(2N  +  l-i>  j-1  (ij  -  Yj)(J -DK21M-1  - J) J 


One  advantage  here,  the  option  of  calculation  in  extended  precision  of  the 
coefficients  of  hence  allowing  one  to  "Monte  Carlo"  the  direct  solu¬ 

tion  and  examine  questions  of  behaviour. 

The  linear  set  of  equations,  previous.  Is  known  ill-conditioned.  In¬ 
tuitively,  the  compressed  visible  range,  relative  to  the  range  of  hopeful 
extension  (say  1 8 1  <  ««  h|b|,  N  -  5,7)  and  the  slow  convergence  of  the 
sampling  series  tend  to  require  a  conditioning  problem.  Ck»e  can't  enlarge 
the  visible  range.  Since  we  can  sample  freely  however  in  its  continuum,  it 
is  reasonable  to  aak  if  the  "sampling"  series  rate  of  convergence  ney  be 


enhanced.  To  this  end  we  discuss  briefly  the  convergence  properties  of  the 
sampling  series  and  produce  the  long  awaited  "quick"  proxlnates  promised 
previously.  He  then  use  these  new  series  to  obtain  an  analytic  continuation 
of  the  power  pattern  function. 


l.e.,  8^  ■  Ai  +  8^ 


IB 


To  start,  consider  the  Man-square  error  in  truncation  of  the  sampl¬ 
ing  series  over  11 .  Recall  the  orthogonality  relationship  of 
(sinc(u-n)  on  R  ,  and  let 

«M(u)  -  f(u)  -  f»(u)  ■  ,  X)  *  f (n)sinc(u  -n) 

M£n+1 

El  * 

—  f (n)sinc(u -n)  and  f  is  of  exponential  type  ir, 
*  |n|<  N  * 

band limited  [-«,«]  .  The  m.s.e. 


/OD  -CD 

CjjTjjdu  »  J  Y*  f (n)sinc(u -n)£  ?(n*)sinc(u -n’)du. 
^  IT  H  11 


Note  orthogonality  kills  the  croes  terms  hence  E  »  — *■  L  | f (n) | 

f  (n)  *  |n|>N  +  l 

The  are  Fourier  coefficients  and  decay  at  least  as 

n 

0(—  -1.  Thus  the  m.s.e.  in  truncation  is  (cf.  [10],  p.  170)  with  some  work 
n 

'v*  .  As  mean-square  error  is  quadratic,  one  might  expect  and  would  find 

N  . 

absolute  error  *  »  hence  the  "tardy"  convergence  observed.  Indeed,  if 


N  »  1 


‘i  •  lf  -  f« 


r  f (n)  sinirt  i  r  If (n)  1  ,  isirnrt 

1  >  N  *  I  -  |nJ  >  n  |t  -  n|  ‘  w2 


£  W -  °<|K  '  /B 


N  -  t 


*  e  . 


*  .  sin  wx 

sine  x  *  ~  •— 
wx 


’for  convenience  2a  »  w,  renormalisation. 
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Published  results  Indicate  [11,  12,  13]  that  though  It  Is  possible  to 
"fine  tune"  the  above,  the  aagnltudes  must  stand*  We  now  derive  saapling 


proxlmates  that  display  enhanced  polntvlse  convergence  of  arbitrary  order, 
given  knowledge  of  f  and  Its  derivatives  in  the  visible  region,  at  arbi¬ 
trary  points  Interior  to  the  sample  points  (unknowns  ,  P(^j)  ). 

Let  g(x)  be  of  exponential  type  v  (this  scaling  to  slapllfy  nota¬ 


tion)  so  that 


g<*)  -  f 


e**%(t)dt 


where  $(t)  is  lr  [-n,n].  Porn  the  quantity  g(s) 


.  g(»)  ~  8(0)  g*(o)g 


where  without  loss  of  generality  we  require  knowledge  of  g,  g*,  g"  (later  at 

z  *  0.  One.  can  readily  generalise  to  any  point  In  the  visible  region.  Mote 

g(z)  has  a  removable  singularity  at  zero  and 

,  ,2 

»<0)  +^«*(0)  4^-  g"(0)  +  •••  -  g(0)  -  *'(0)^* 

11m  g(z)  -  11m  - - - * — 


z  0 


z  -*•  0 


These  quantities  exist  and  will  be  Ir  (R)  for  g(z)  of  modest  growth  i.e. 
g(z)  t  0(z°)  where  o  <  3/2.  In  this  case,  g(a)  l^dl)*  so  by  the  Paley- 
Wlener  Theorem*,  since  g  Is  of  exponential  type  w,  it  is  bend -limited 
[-*,*].  Now  replace  g(z)  by  Its  Cardinal  (Sampling)  Series  representation* 


*<*>  •  £  8inwtl  “  *2^ "  V21  ]* 


We  require  the  following  identities; 


nw(z  -n) 


1  sin* 


-  i  (1  -  ~~)  where  hereafter 


*  2 
The  entire^ function  f(z)  is  of  order  1,  type  e  and  belongs  to  1*  (R) 

iff  f (s)  -  j  eic**F(s)ds  where  F(s)  is  l.2 (-1,1)  . 
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the  prime  will  denote  deletion  of  n  ■  0  term.) 


-  li*  7U-I  -  lim  4^-0, 

i  4  0  1  ■5I 

Is  bounded, 

2 

1  (®)>  end  of  exponential  type  it,  we  have,  again  by  the  Paley- 
Vttener  Th,,  the  entire  quantity  is  band-limited.  Thus  using  again  the 
sampling  theorem 


Proof:  Observe  lim  i  (1  - 

*  -  0  *  ” 

hence  the  quantity  is  entire.  As  sinc(z) 


Avi  sinirts  _  V'  A  tt  *iniT\  (ain»rz)(-)n  ..  r"  1  (-)nsinirz 
zU  its  >  ^  n  U  -  *z  >  i(I50 - ^  * 

•  5(0, n) 

n  2 

«-«  1  sinirz(-)  1  ,,  sinffz,  1  linn 

Ui;  2  v(*-n)  "  2  U  '  its  1  "  3!  n  “ 

n  z 


Proof:  Consider  [1  -  1,  note  the  quantity  has  a  removable 

2 

singularity  at  zero  and  is  thus  Entire,  of  type  *.  Further,  the 

2  2 

limiting  value  at  zero  is  ir  /3!  *  As  the  quantity  is  also  1,  (It)  the 
sampling  theorem  applies 


slrarz 

irz 


w 


All 


sin  nr  .  sin7»(-)n  . 

rn  J  v(z  -n)  T  3! 


ainwz 

wz 


These  results  set  directly  into  the  sampling  expression  for  g(z), 
holding  in  mind  that  g(0)  •  g"(0)/2  give: 

t(»  -  -fjp  ,”(0)  .  2'iiWsjUEl  .  ,.(0)  i(l  .  •gl*1 

W  wn  (z  -  n) 


2 

_  zA\r  1  im  siwrz  1  w  simrz. 
,(o)t-j<i  3i  7.—: 1 
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« 

Now,  writing  g(z)  in  terns  of  g(z)  as  previously  defined  gives 


g(s)  -  g"<0)  +  g’<P)  +  g(0)  [-2^p  +  ] 


+  £  '  C~)°  *  sinirz  »  g(n)  »  z 


» 


n  •  it  •  (z  -n) 


The  result  is  confirmed*  [14]  In  the  literature.  He  note  that  knowledge  of 
higher  derivatives  at  a  point  in  the  visible  region  produces  a  ’’faster" 
family  of  series  [15].  Also,  estimating  the  (second)  derivative  In  the 
presence  of  noise  is  difficult  enough  -  again  an  "ill -posed"  problem,  but 
a  well  studied  one.  (e.g.  splines,  interpolations  etc.  [16]).  To  obtain 
these  even  "faster"  estimators  one  need  only  substitute  more  terms  of  the 
Taylor  series  cf.  g(z).  Late  literature  quotes  a  "slower"  result  requiring 
knowledge  of  only  one  derivative  [17]  -  the  proof  does  not  generalize: 

g(z)  .  *  g’(0)  +  g(Q) sinirz  +  *7*  (-)nsin  w*  »  g(n)  •  t 

*  **  V  *n(z  -  n) 

Somewhat  less  ambitious  but  more  practicable.  At  this  juncture  it  would  be 
possible  to  establish  error  bounds  -  a  little  thought, however,  would  convince 


* 


not  shown,  general  discussion  [15 ]  gives  gent  for  proof. 
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the  reader  that  the  previous  bounds  are  improved  by  a  factor  of  a  in 
the  series  denominators  so  results  are  one  to  two  magnitudes  better. 

We  illustrate  for  S(t)  =  |sint|  the  explicit  series  fora.  (A 
numerical  table  of  convergences  and  error  is  in  Appendix  B,  pt.  II).  Bote 
|sint|  has  the  required  vanishing  properties  at  +  tt  and  is  not  too 
smooth. 

«(«)-  f  «l!t  lilnthn  -  e(0) .  4 

J-r  t  l  -  i 


g(n) 


n  odd 


g(n)  »  V(1  -  n  )  n  even. 

A  simple  way  to  find  derivatives  is  the  direct  differentiation  of  the 
integral  representation.  Whereupon,  g’(0)  =  0,  g"(0)  ■  8  -  2n2  . 
Thus  the  three  series  interpolations  for  g(s) 


1)  g(*)  -  L 


simrz  ,  k  x 

«.»  d  ^  (rr7) 


(Cardinal  Series) 


2)  *b).H£SL*  E 


(eimr*)  •  z(- 


TTZ 


even  n 


1-  n 
t m(z  -  n) 


simre 


1b' 


3)  «(■)-  [8-^Sl£2&+4[5^  +  2!|i2!!£  1+  £'  —JL  =JLi 

^  ™  3-  even  n  n%(z  -  n) 


Bote  the  series  are  0(-rr),  0(-qr),  O(-tt)  convergent  respectively.  However, 

h^  nH  tr 

large  values  of  s  require  a  "cancellation"  convergence  -  the  subtraction 
of  large  quantities  from  each  other  -  which  can  get  sticky*  We  will  be 
operating,  however,  with  z  visible,  so  given  accurate  measurement  this 


fact  is  of  little  import. 

Heedless  to  say,  one  can  now  rework  the  solution  of  the  truncated 
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linear  sets  previously  presented  using  these  "faster"  sampling  series.  (A 
typical  example  is  discussed  in  Appendix  II).  The  analytic  arguments  are 
paralleled  by  those  previously  given  and  we  omit  the  details.  It  is  easy 
to  see  that  the  double-alternant  form  is  preserved  in  the  "faster"  equation 
set  and,  it  is  possible  to  consider  the  behaviour  of  the  exact  solution. 

In  passing  we  note  that  these  "faster"  formulae  for  band-limited  functions 
may  be  used  in  the  previous  integral  expression  for  the  mutual  coupling  to 
derive  a  rapidly  convergent  expression  for  the  coupling  involving  only  the 
sample  points.  Consider  the  expression  for  g(z)  [— tt, tt]  with 

g'(0)  =  0  (maximum  at  0°  «  8)  as  an  example.  (We  use  the  second  expres¬ 
sion  requiring  knowledge  of  the  first  derivative) 

rg(g)slmru  +  y,  (-)n(sinyu)u-g(n),c-luDdu 
iuD  f»>  ttu  7m(u  -  n)  3 


%  z  .  r 

CP  21  -  I 


g(u)e" 


J? T7 


du 


/ 

•T-ii 


JFT? 


The  first  integral  has  been  evaluated  as  a  convolution.  The  remaining  may 
be  easily  found  with  the  following 


4gSiB*£l 

TD 


f  °°  .  .  .  -iuD  .  simr(u-n) 
jL  ttu  -  mr 


du 


which  is  directly  proportional  to  the  required  quantity.  (Note  powers: 

u11  «— »&^n\x)  etc.).  The  numerical  evaluation  of  3nI  is  simple,  Just  the 

d 

introduction  of  a  cot  6  term  in  the  integrand  of  the  1~  form  of 
I(n,B,D)  (cf.  p.  7)»  So 

slnc(u-n)eiuD*  udu  •  (&ffi-) 


-  g(o)«o,,,D)  ♦  E'  1  [|$niPiMi 


JT- 


2k 


where  the  series  convergence  is  now  —  versus  —k  previously.  That 

1  “ 

SD*D  ~  n  is  obvious  fran  the  convolution  form  of  l(n,  g,D)  over  a  finite 

interval,  i.e.the  Fourier  coefficients  of  the  's  first  derivative. 

These  arguments,  again,  may  be  made  quite  general  with  the  accurate 
knowledge  of  each  successive  derivative  of  g(u)  at  a  point  enhancing  the 
rate  of  convergence  of  Z21  by  an  order  of  magnitude  as  above.  Thus  the 
series  form  of  is  made  practicable  if  it  is  possible  to  exprapolate 

the  first  few  values  of  P(^)  outside  the  visible  region.  A  perusal  of 
the  appendices  hints  that  while  least  square  "ball-park"  results  are 
readily  forthcoming  much  work  needs  yet  be  done.  The  stabilization  of  an 
ill-posed  problem,  described  heretofore  in  totally  analytic  terns,  e.g. 
simultaneous  solution  of  linear  equations,  may  be  accomplished  if  constraints 
exterior  to  the  data  are  introduced.  Now  we  leave  the  discussion  of  the 
" sampling -proximates"  solutions  to  explore  the  general  methods  of  Miller  - 
Tikhonov  regularization  for  such  problems.  We  will  later  return  to 
these  truncated  linear  sets  as  subject  for  regularization. 

Regularization  is  necessary  when  the  problem  at  hand  is  "ill-posed" 
which  is  not  to  say  unphysical.  Required  in  this  context  may  be  1)  numeri¬ 
cal  differentiation  2)  analytic  continuation  3)  solution  of  Fredholm  or 
Volterra  equations  of  the  first  kind  etc.  Many  proper  physical  situations, 
e.g. computation  of  atmospheric  density  profiles,  nuclear  scattering,  inverse 
scattering,  remote  sensing,  and  image  enhancements  in  optics  give  rise 
to  models  that  are  "ill-posed"  or  "conditioned".  The  problem  at  hand, 
numerical  analytic  continuation  of  bandllmited  functions,  may  be  seen  as 
ill-conditioned  with  the  aid  of  the  Riemann-Lebesgue  Leona.  In  the  Fredholm 
formulation,  we  wish  to  solve  for  f(x),  with  g(x)  *  /  f(y)elyxdy 

4  vx 

g(x)  for  x  c  v  given.  Nov  e  is  a  smooth  kernel,  thus  this  integral 
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operator  does  not  have  a  bounded  Inverse.  To  see  this,  let 

f  (y)  =  sin  n  y  ,  note  that  for  n  big  |  /  f  (y)e  may  be  Bade 

€  €  €  € 

as  small  as  desired  by  the  Riemann-Lebesgue  Leona.  Hence  f(y)  and 
f(y)  +  f£y)  are  indistinguishable  as  "solutions"  i.e.  finite  changes  in 
f(y)  may  produce  infinitesimal  changes  in  g(x).  However,  if  one  restricts 
suitably  the  characteristics  of  the  solution,  f(y),  e.g.  take 
/  (f ' (y))2dy  a  minimus,  one  produces  a  restricted  problem  whose  solu- 
tion  is  unique,  with  the  continuous  dependence  on  data,  g(x),  restored. 

To  this  end  we  consider  three  such  formulations.  It  is  not  known  as  yet 
which  of  these  procedures  are  optimal  or  feasible.  This  must  be  the  sub¬ 
ject  of  future  msnerical  research.  We  can  only  outline  the  various  methods 
and  indicate  our  reasons  for  preference  or  otherwise. 

Tikhonov  regularization  [16,  I9,  20]  may  be  readily  applied  to  the 
collocative  solution  of  the  Fredholm  equation  of  the  first  kind  [21]. 


In  such  a  solution,  the  required  extrapolation  is  Just  a  matter  of  quad, 
rature  i.e.  given  g(t)  =  f(x)e1,*xdx  and  given  f(x),  g(z)  is 
"defined"  by  quadrature  outside  V  =  [-8,0].  The  problem  is  the  determina¬ 
tion  of  f(x).  Let  Tf(x)  *  f  f(x)eizxdx  and  take  an  12(V)  setting 
[21].  If  g  is  the  exact  pattern  and  gy  the  measured  pattern,  let 
Iky  -  8l|y  < v  .  Introduce  the  functional  Mffjg^or)  =  llTf-gJ2  +«W(f) 

where  W(f)  *  C0||f||2  +  cillfil|2  and  ||b(x)||2  *  J h(x)E(x)dx  ,  i.e.  the 

2  ^ 
standard  1.  (V)  norm,  or,  the  regularization  parameter,  is  >  0,  with 

CQ, likewise,  A  theorem  [18,  Chap.  II]  of  Tikhonov  requires  that  for 

o 

all  gy  «  1  (V)  there  exist  a  unique,  continuous,  differentiable  foy 
minimising  M  given  or,  gy.  This  f^y  is  seen  to  be  a  stable  solution 
to  the  original  problem  in  least  squares  sense  with  f^  -» f  as 
{»,v}  -*  0  [18]  •  The  orW(f)  acts  as  a  sort  of  Lagrange  multiplier,  restricting 
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oscillatory  or  excessive  behaviour  of  f ,  frequently  or  is  quite  . 
providing  "Just  a  bit"  of  damping.  For  such  an  f^,  the  first  variation 
must  vanish  i.e. 


where  T*f 


6f[M(f  ,8^001  =  0  4= 
=  f  e_izxf(x)dx  and 

•'wir 


[T*T  +  *(C0I  -  ^D2)]^  =  T*gv 


™  -fv  ««•>*»' 

p 

Note  T*  is  the  adjoint  on  ]L  (V) ,  D  the  differentiation  operator,  I 
identity  etc.  Thus  a  formal  solution  is 

f  =  L'Vgy  where  L  =  [T*T  +  <r(CoI  -  C^D2)  ] 


Note  that  the  choice  of  or  has  not  been  specified.  Here  seme  of  the 
"ad  hoc"  flavour  of  this  technique  sets  in.  For  specific  values  of 
Cq  =  1,  =  0  [21]  it  can  be  shown  that  IJf^  -  fj|  <  v/S^r  where 

fo  minimizes  M(f,g,a)  (noise-free  smoothing  .  For  {|  f  a  —  f  ||<y  (|f|| 

where  is  max  |>^|  for  T*Tf  «  "K^t  with  y «  ||f||,  a  <  7^||f  ||.  Hence 

with  some  knowledge  of  ||f |j^ J | current | Magnitude  in  "slot")  one  has  an 
upper  bound  on  or.  A  bound  on  Ijf^  -f ||  may  be  calculated  and  [21,19] 
given,  thus  ([f^  -  f||  <  v/2 Jx  +  (a/a+  >^)||fj|.  If  the  self -Impedance  of  the 
antenna  producing  the  given  pattern  is  known,  one  might  pick  a  such  that 
the  pattern  solution  provides  the  best  approximation  to  the  measured  self- 
impedance  (to  be  evaluated  by  the  integral  forms  involving  visible  and  in¬ 
visible  patterns  for  self -Impedance  discussed  previously.)  In  general,  the 
selection  of  or  is  a  "questionable"  procedure  subject  to  nvaerical  experi¬ 
ment.  In  using  the  Fredholm  approach  to  continuation,  we  are  required  to 


solve  for  a  quantity,  f,  which  is  then  integrated  with  another  kernel, 
providing  the  required  extrapolation.  One  might  expect  that  the  reinte¬ 
gration  would  obliterate  some  of  the  error  in  solution.  However,  we  are 
also  finding  the  inverse  of  a  transform  ,  f,  only  to  retransform  it,  a 
seemingly  redundant  procedure.  In  fact,  it  can  be  shown  [22,23]  that  the 
error  In  finding  f  is  subject  to  a  logarithmic  continuity,  whereas  the 
optimal  extrapolation  error  for  extracting  g(u),  u  c  invisible,  is 
subject  to  a  Holder  continuity  i.e.  error  in  |f  -  f  |  ~  (An  v)~^  as 

*v  VIsl Ue»  g  Whereas  l6^)-8v^)h  (v)®  as  gy g  with 

u  €  invisible  region 

0  <  p  <  1  for  optimal  methods  in  both  regards.  The  hazards  of  logarithmic 
continuity  need  numerical  investigation  in  this  connection.  An  important 
benefit  in  the  use  of  the  Tikhonov  scheme  is  its  ease  of  Implementation. 
Whereas  "optimal"  methods  referred  to  above  make  use  of  the  Prolate  Spher¬ 
oidal  Harmonics,  the  procedure  outlined  here  requires  only  simple  quadra¬ 
tures  and  the  equations  in  question  are  stable.  The  descreticatian  of  all 
operators  with  any  Gaussian  Quadrature  scheme  (to  minimize  required 
sampling  of  the  visible  pattern)  presents  no  nuisance  and  the  inversion  of 
L  presents  no  difficulty.  (Also  numerically  cheap: )  Details  of  numerical, 
experiments  [21,16]  are  also  encouraging. 

A  more  elaborate  technique  of  regularization,  due  to  Miller  [2k]  and 
extended  by  Viano  et  al.,  [23,  25  see  also  the  exhaustive  references  listed 
there]  makes  good  use  of  the  Prolate  Spheroidal  Harmonics  to  establish  the 
analytic  continuation  of  bandllmited  functions.  While  quite  elegant,  this 
theory’s  use  of  P.S.H.  condemns  one,  in  trying  to  implement  it,  to  a  host 
of  numerical  nightmares.  Even  recent  algorithms  [261  for  computation  of 
P.S.H.  eigenfunctions  and  eigenvalues  require  messy  matrix  bisection/Stnmi 
e.v.  sequences  procedures  which  allow  no  simple  chax^e  or  rescaling  without 
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complete  reccmputation  at  great  expense.  The  only  virtue  in  these  techni¬ 
ques  is  their  "best -possible"  nature.  Some  consideration  shovs  the  method 
instinctively  no  different  from  those  suggested  [1]  in  our  first  progress 
report.  We  consider  briefly  these  results. 

Let  X  denote  the  set  of  all  permissible  pattern  functions  (i.e. 
f  e  X  is  band-limited  appropriately,  /  (f  |2du  <  «°  etc).  For  f  e  X 
let  Af  denote  whatsoever  is  subject  to  measurement  (the  "perfect"  data 
in  V)  approximately.  Take  h  as  the  measured  data.  Af  and  h  are 

O 

contained  in  the  "data  space"  Y  (in  this  case  1  (V)  ).  To  stabilize 

the  problem,  one  constrains  f  by  "boundary"  or  "global"  values.  Let 

B  be  a  linear  operator.  Bf  provides  some  constraint  operation  on  f. 

Bf  e  Z  the  "constraint"  space.  Then  if  f  satisfies 

o 

II  Afo  —  till  Y  —  *  (sat*8*1-68  expected  data  fit) 

||  BfQ  ||  z  £  E  (satisfies  global  bound) 

it  may  be  shown  fQ  also  satisfies  (by  Lagrangian  arguments) 

II  AfQ  —  h II  y  +  ^  II  Bf  ||  ^  <  2f  This  last  equation  may  be  mini¬ 

mized  by  solution  of  the  normal  equations  to  obtain  an  f  which  can  be 
shown  to  be  "optimally"  fit  to  fQ  independent  of  norm  used  to  measure  the 
error  [see,  e.g.,  24],  If  the  criterion  of  error  is  <  •  >g  ,  it  may  be 
shown  <f  -  fQ>e  <  ^2  M(*,E),  where  M(«,E)  is  defined  as 
sup(<x>e  |x  e  X,  ||  Ax ||  y  <  «  ,  ||  Bxll^  <  E}  -  the  "best  possible  stability 
estimate."  Called  as  such  since  M(  e,E)  gives  the  "size"  of  all  f  that 
aatisfy  the  global  and  data  bounds,  for  such  (f^^),  ||A(f ^  -  fj)  II y  <  2  , 

(^(fj  -  f 2)  ||  <_  2E  so  that  <fj  -  f2>e  <_  2M(  e  ,E)  [see  23,24)  etc. 

Whence  we  may  take  f  -  [A*A  +  (fr)^  B*B]”*A*h  as  a  regularised 

C 

solution  whose  error  (in  <  >@)  is  "best  possible"  but  for  a  factor  of 
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^2.  Numerous  estimates  may  be  fashioned  in  terms  of  the  eigenfunctions 
and  spectrum  of  A,  A*A,  etc.  for  the  various  error  norms  [cf.  23]  but 
these  estimates  generally  require  full  knowledge  of  the  spectral  de¬ 
composition  of  the  various  operators  -  something  not  readily  available 
here. 

Our  problem  is  placed  in  the  above  setting  when  one  considers  the 
problem  of  extrapolation  of  optical  images  (see  [25]  for  physical  details.) 
We  require  explicity  that  the  object  image  Is  of  bounded  energy  E  in 
TL2  [-1,1]. 

2 

Let  x  »  x(t),  t  €  K  ,  be  L  (R)  ,  bandlimited  such  that  the 

image,  x(t)  »  I  c^v  e2iritu  an(j  take  Ax  to  be  the  restriction 

J- c/2nir 

of  x  to  the  interval  [-1,1].  (We  assume  that  the  Images  are  known,  in 
a  system  of  unit  magnification,  over  an  interval  equal  to  the  support  of 
the  object,  [-1,1].)  Take  h  -  Ax  +  z,  the  image  corrupted  by  measure¬ 
ment  error.  The  problem  to  estimate  x(t)  given  h(t).  In  this  context 
we  identify  the  following 

1)  X 

2)  Y 

3)  Ax(t)  «  x(t)  t  €  [-1,1],  ■  0  otherwise 

a)  A*y(t)  -  J  y(s>ds  y  eY  etc. 


« - *  {H*( 


delimited  f} 


functions  of  compact  support  [-1,1]  ^  1L  (-1^ 


It  is  possible  to  specify  a  stabilization  constraint 


x(t)  ■  j  ~S^v(s)ds  where  J  |v(s)|^ds  <  Z2 

(band-pass  filter  etc.) 
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with  regards  to  optical  object  energy.  With  suitable  normalization 
one  can  take  B  ■  (A*)”*  and  after  much  dust  get  a  "best"  approxi¬ 
mation 


x<t)  “  £ 


k“0  >k+<t>2 


Vk<c»t) 


i  r1 

with  hfc  «  (h.v^Y  *  — —  J  h(t)^(c,t)dt  and  ^  the  appropriate 
P.S.H.  This  review  is  not  designed  to  be  complete  -  details  may  be 
found  in  [23,24],  Others  have  considered  the  use  of  P.S.H.  to 
extrapolate  [27]  band-limited  functions  with  constraints  and  achieved 
similar  results.  Though  keenly  analytic  in  nature,  they  all  evolve  as 
quite  cumbersome  due  to  the  use  of  P.S.H.  We  therefore  continue  to 
ignore  such  schemes  save  as  possibly  benchmark  measures  for  other  pro¬ 
posals.  It  shoud  be  noted  that  the  constraining  operator  B,  data  fit 
constants  (c,  E)  need  not  both  be  known  [24]  to  use  the  Miller  method. 
Through  certain  concavity  properties  it  suffices  to  know  just  one  of 
these  quantities.  The  selection  of  B,  the  linear  constraint  operator,  is 
quite  open.  Further  study,  elsewhere,  would  seem  appropriate. 

Let  us  note  that  knowledge  of  the  self  impedance  of  a  given  antenna 

does  not  qualify  as  a  constraint  in  the  sense  given  above.  (This  does 

not  exclude  its  use  as  a  Lagrangian  constraint  as  we  shall  see  -  it  simply 

will  not  allow  us  to  use  it  to  stabilize  the  ill -posed  problem.) 

Specifically,  the  input  impedance  Z.,  a  /  — i8  a  functional 

ln  ,  /-»  x  ”  - 


of  P(u)  unless  we  are  allowed  to  vary  B.  “  But  P(u)  implicitly 

d 

(especially  in  3  )  contains  £,  hence  rendering  different  frequency  measure¬ 
ments,  unless  P(u)  is  B  invariable,  worthless.  The  question  of  what 
information  may  be  had  from  different  freqyebct  measurements  of  the 
visible  pattern  and  impedance  is  open.  It  would  seem  necessary  to  have, 
a  priori,  some  notion  of  antenna  geometry  to  pursue  this 


'2  2 
6  -  u 
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problem. 


That  self-impedance  specification  will  not  stabilize  the  prco  -a  is 
2 

seen  In  the  previous  1L  (K)  context  of  the  Miller  example.  That 
1 


<  P(u),- 


—  2  /<|p\  A2  -  u" 

1L  (]R)  vector 


>  »  Const  merely  implies  the  projection  of  P(u)  on  the 


is  fixed  -  this  does  not  preclude  wild 


oscillations  as  previously  discussed  over  sets  of  finite  measure  as  would, 

say  some  limit  on  the  average  value  of  P'(u)  (mean  squared).  This  does 

not  mean  one  ought  discard  this  information  but,  that  it  is  not  enough 

hence,  one  can  use  it  as  another  Lagranglan  constraint.  (It  turns  out  that 

if  one  tries  to  continue  the  pattern  function  (not  power)  the  Impedance 

f  If  I  ^ 

operator  can  be  reformulated  to  be  of  service,  /  — -*■  ■* —  du  places 

Jn 


some  constraints  on 


L 


Ifl^du-  but  this  case  will  not  admit  useful 


mea surement.) 

The  final  method  of  regularization  applies  to  the  solution  of  linear 
equations  with  error  in  the  "right  hand  side"  (i.e.  the  specific  measure¬ 
ments  of  g(^)  +  Error  (n))  and  the  stabilization  thereof,  [28,29,18, 
Chap.  Ill],  These  results  apply  directly  to  the  truncated  linear  sets 
considered  in  the  text.  The  setup  is  as  follows: 


given  Af  ■  g  +  e  and  £  <  e^  «  1 

how  can  one  vary  €  (since  A,  g  are  fixed  this  is  equivalent  to 

"vary  f")to  minimize  the  a  quadratic  constraint  on  |, 

Q  *  X.  h4.f .f.  •  f*Hf.  There  are  two  cases  of  interest: 

ij  1J  1  J  - 

1)  there  exist  {g^}^  such  that  Q  is  minimized  and 

II  H  ”  8 if  <  •  In  this  case  g^  is  a  variational  solution 

to  6Q  ■  0  and  Af  -  g  +  *  -  which  is  unlikely  - 

5 

i  »  index  of  vectors. 
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is  however  a  valid  solution. 


2)  For  Q  such  that  a  variational  minimum  does  not  occur  within 

2  2 

{f}  such  that  Hell  <_  e  ,  we  know  that  the  minimum  must  occur 
on  the  boundary  of  {f  | (Af  ■  g  +  e  or  f  ■  A  *(g  +  «))  and 
||  l||  -  e  }  or  simply  for  {f  |  ||e  ||  -  e). 

Therefore  we  may  use  the  method  of  Lagrange  Multiplier(s)  to  solve  the 

following  equivalent  problem: 

Minimize  Q  ■  f*Hf,  subject  to  the  constraint 
-  «- 

(Af  -  g)*(Af  -  g)  ■  e^  i.e.  minimize  the  functional 

f  *Hf  +  y  €  *  t  . 

•  *-  o  “  - 

In  our  context,  we  may  require  a  second  Lagrange  Multiplier  to  take  into 
account  any  knowledge  of  the  self -impedance  associated  with  the  measured 
visible  pattern.  In  general,  one  may  also  require  “  f  •  c. 

Thus  for  F(f )  ■  f*Hf  +  y  +  y,  f*c  we  require  6,F(f)  »  0 

i.e.  {[f*A*Af  -  g*Af  -  f *A*g  +  g*gW  +  f*Hf  +  Yif*c)  -  0 

3fk  -  .  2  -  -  s  5  .  -  o  ...  1-  - 

or  (2A*Af  -  A*g  -  A*g)y  +  2Hf  +  y.c  ■  0  i.e. 

as*  s  2  -So  s-  I-  - 

(**^0  +  S>i  *  Yof?  -  T fi 

where  y  y.  are  to  be  selected  to  satisfy  the  discrepancy  and  impedance 
o  1 

constraints,  (In  practice  the  second  constraint  is  linear  and  it  incor¬ 
porates  directly  into  the  solution  of  the  above  equations;  if  A,  f  are 

S»  ** 

n  x  n,  n  etc.  ,  we  have,  with  the  equation  Z^  ■  f  •  c  n  +  1  linear 
equations  in  terms  of  the  n  +  1  unknowns  (f  and  y^  •  using  yQ 


« 


33 


2 

'hs  a  known"  parameter  which  we  vary  to  (experimentally) satisfy  **€■  e  ). 
The  application  of  this  technique  to  the  previous  sampling  equation  set(s) 
is  immediate  and  will  be  considered  and  coded  as  soon  as  time  permits. 

A  few  remarks  concerning  the  method  are  in  order. 

1)  The  quadratic  constraining  function  Q  is  essentially  arbitrary. 

Possible  forms  for  H  are  "minimum  variance",  first  difference  (tabular), 

s 

and  prior  departure  from  given  form,  e.g.  the  asymptotic  pattern  form.* 
Which  form  or  forms  are  suitable  must,  perforce,  be  found  in  an  empirical 
fashion. 


*  2)  One  can  show  that  the  variation  of  f  with  y  is  slow  and 

W  «  O 

2 

the  ensuing  numerical  satisfaction  of  f  *  f  *  e*  is  "easy"  with  convex¬ 
ity  arguments. 

3)  The  method  is  a  "generalized"  least  squares  technique  no  more 
difficult  than  the  solution  of  the  normal  equations.  That  it  is  possible 
to  satisfy  and  solve  the  Lagranglan  constraints  (no  Inconsistency)  is 
not  clear  from  the  above  arguments  -  we  refer  to  [16]  for  a  discussion 
of  this.  The  method  was  one  of  the  first  tractable  solutions  to  the 
numerical  problem  of  the  First  Kind  Fredholm  Equation,  [29,  circa  1962]. 
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S 
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Conclusions  and  Extensions. 


We  have  presented  the  rudiments  of  a  theory  for  mutual  coupling 
between  canonical  minimum  scattering  antennae  requiring  knowledge  of  the 
total  pattern  function.  The  problem  of  extending  the  visible  pattern 
function  throughout  the  entire  required  domain  was  solved  by  introducing 
the  constraint  of  band-limitation.  Whereupon,  a  series  form  representation 


of  the  required  coupling  integrals,  utilizing  various  extensions  of  the 


*note  g(f)  ■  f  e**xf(x)dx  * 

J-n 

an  "observable  "  quantity. 


see,  e.g.  [30], 
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sampling  theorem,  was  developed.  Lastly,  we  have  proposed  several 
techniques  that  should  allow  the  computation  of  the  invisible  pattern 
values  from  the  available  data.  Several  well  known  regularization 
techniques  were  drawn  from  the  literature  in  the  context  of  our  pro¬ 
blem.  The  present  research  must  now  shift  emphasis  to  the  numerical 
validation  of  the  "sampling"  forms  of  continuation  and/or  comparison  to 
the  Miller -Tikhonov  method.  The  ill  conditioned  numerical  problem  at 
hand  required  an  extensive  investigation  into  possible  analytic  forms  of 
solution  and  having  accomplished  this,  we  now  hope  to  devote  substantially 
more  time  to  the  actual  computations  required.  The  appendix  confirms 
at  least  order  of  magnitude  results  are  possible  so,  at  worst,  "reasonable" 
bounds  on  the  couplings  may  be  expected. 

There  are  several  points  novel  to  our  approach.  Intrinsic  is  the 
restriction  to  coupling  amongst  like  antennae.  This  is  required  by  the 
difficulty  of  phase  measurement.  The  problem  of  phase  retrieval  Ie.g.  31] 
is  well  known  in  optics  but  usually  requires,  say  for  a  Logarithmic 
Hilbert  transform  resolution,  knowledge  of  the  pattern  magnitude  and  its 
zeros  over  a  relatively  large  interval, information  here  unknown.  It  may 
be  possible  to  model  or  bound  phase  effects  on  the  coupling  and  extrapola¬ 
tion  formulae  previously  obtained,  whereby  the  generality  of  the  analysis 
would  be  restored.  However,  save  at  the  expense  of  accurate  phase  measure¬ 
ment,  the  problem  remains.  The  most  useful  observation  is  the  bandllmlted 
nature  of  the  pattern  power  function.  This  lone  fact  restores  the  accu¬ 
racy  needed  to  consider  extrapolation  seriously*  With  the  hope  of  extend¬ 
ing  the  techniques  to  multiple  dimensional  antennae,  we  have  kept  the 
sampling  theorem^  multidimensional  character  In  mind  [32].  These  notions 
would  hopefully  allow  analysis  of  measured  couplings  between  planar  aper¬ 
tures  and  the  like  to  follow  through,  whereas  the  various  integral  form- 
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latlons  (Fredholm)  would  not  readily  reformulate.  Note  the  assumption  of 

bandlimitation  to  a  degree  holds  Implicitly  some  knowledge  of  that  degree. 

To  the  extent  all  physical  systems  behold  finite  measure  of  energy  and 

d 

geometry  (explicit  in  2~  )  suggestive  of  the  bandllmit,  one  may  be  for¬ 
given  its  accurate  assumption.  The  sensitivity  of  this  assumption  to 
error  is  another  matter.  Lastly,  the  optics  literature,  in  considering 
the  problems  of  image-processing  and  restoration  -  likewise  demanding 
extrapolation  -  may  yet  be  of  further  service. 
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Appendix  A 


Numerical  Illustration  of  the  First  Variational  Method. 


Correct  end  generated  solutions  to  the  extrapolation  problem  for 

various  measure  ere  presented.  For  this  and  subsequent  examples  the 

d 

following  band limited  power  pattern,  typical  of  a  T~  slot  radiator, 
will  be  employed.  Note  the  pattern  is  real  and  even,  whereby  many 
simplifications  ensue. 
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Taking  an  1L  (-n,ir)  setting,  pick  as  the  unknown  proximate  to  be  matched 

on  (-8,8)  A(u)  -  aft  and  invoke  the  constraint  that  aR  ■  a^, 
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where  aQ  -  1  -  P(0),  which  is  known  a  priori.  As  previously  discussed,  the 
{a^}  need  minimise  the  functional 
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The  last  constraints  provide  the  first  two  of  H  +  2  linear  equations  for  the 
unknowns  (a^,),  Y^»  Y2 
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Appendix  B.  Pt.  I 


The  solution  of  Least -Square  Sampling  Proxlaates.  Slow  and  Fa at. 

For  various  values  of  (6,N)  a  least  square  and  direct  solution  of 
the  Harris-type  extrapolation  procedure  la  provided.  Again,  the  even 
slot-llke  pattern  Is  used.  Note  the  order  of  magnitude  "contour"  agree¬ 
ment. 
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(Sampling  series  proximate) 


2N  2 

z_x  v*  stnir*  4*  .  tt  2izslnirc  .  .  ,sln*z  .  zvsinzz, 

*2(l)  "  h  T~ - ~ — TT  +  14  *  *  - +  4l“iT“  +  — ^ - 1 


-2N  n2w(*-n)(l-n2) 


3! 


(2nd  derivative  proximate) 


evaluated  for  2N  •  2,  12,  24 #  36,  48 

at  pointa  z  -  14.0  +  (.221333)*k  ,  k  -  1,2, 3,4, 5 

Note  the  improvement  possible  for  even  g(z)  as  follows! 

|  -£Sr .« ■  .(o)  ♦  L  .» ♦  | 

M 

■  g(0)  +  V  ( — — - )  converging  now  0(*i) 

n-1  *  2  2  n 

i  “  n 


It  la  possible  to  reduce  our  labor  considerably  for  g(n) 


even  or  odd  thus 


ANSWER 


-0.1756904 E-01 
-0.1136048E-01 
-0.47402 68E-02 
-0. 582084 5E-03 
-0.48  95511E-03 

0. 18347 01E-01  0.144414 9E  00 

-0.168 5422E-02  0.8106232E-04 

-0.17  98547E-01  -0.1764774E-01 

-0.176892 9E-01  -0.1755905E-01 

-0.17  61997E-01  -0.1754  7  61E-01 

NEXT  K 

0.2778002E-01  0.2243681E  00 

-0.1989841E-02  0.7362366E-03 

-0.1201773E-01  -0.11 51848E-01 

-0.1154983E-01  -0.1137543E-01 

-0.1144117E-01  -0. 113554 0E-01 

KEXT  K 

0. 242307 5E-01  0.2007780E  00 

-0.1306236E-02  0.1142502E-02 

-0. 533711 9E-02  -0.4887 581E-02 

-0.4911125E-02  -0.4752159E-02 

-0.4812896E-02  -0.4734039E-02 

NEXT  K 

0.9681560E-02  0.8227 921E-01 

-0.369913 9E-03  0.6351471E-03 

-0.83 032 62E-03  -0.647 544 9E-03 

-0. 652607 5E-03  -0.5903244E-03 

-0. 6120130E-03  -0.5817413E-03 

NEXT  K 

-0.8909047E-02  -0.7764 53 0E-01 

0.21534 04 E-03  -0.7429123E-03 

-0.2517626E-03  -0.4310608E-03 

-0.4224665E-03  -0.4873276E-03 

-0.46111 64E-03  -0.4949570E-03 


